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We investigate the correlations in the current fluctuations 
at different terminals of metallic diffusive conductors. We 
start from scattering matrix expressions for the shot noise 
and use the Fisher-Lee relation in combination with diagram 
technique to evaluate the noise correlations. Of particular 
interest are exchange (interference) effects analogous to the 
Hanbury Brown-Twiss effect in optics. We find that the ex- 
change effect exists in the ensemble averaged current corre- 
lations. Depending on the geometry, it might have the same 
magnitude as the mean square current fiuctuations of the shot 
noise. The approach which we use is first applied to present a 
novel derivation of the 1/3-suppression of shot noise in a two- 
terminal geometry, valid for an arbitrary relation between the 
length and wire width. We find that in all geometries corre- 
lations are insensitive to dephasing. 

PACS numbers: 72.70-Hm,73.23.Ps,73.50.Td 



I. INTRODUCTION 

The shot noise in mesoscopic systems |1[ continues to 
attract the attention of both theorists and experimen- 
talists. For diffusive conductors, which are considered 
here, the two-terminal shot noise is quite well studied. 
The spectacular 1/3-suppression of the shot noise with 
respect to the Poisson value, 

S{lj = 0) = r,eGV 

(here, as usual, S{llj) is the Fourier transform of the 
current-current correlator, S{t) — (A/(i)A/(0)), while 
G and V are the conductance of the wire and the applied 
voltage, respectively; A/ = I{t) — (/)), was derived in 
three different ways: from the distribution of transmis- 
sion eigenvalues in a wire H], semi-classically from the 
Langevin equation H, and through a microscopic cal- 
culation of local current densities [0. Later, Nazarov 
g claimed that this 1/3 suppression holds for an arbi- 
trary two-terminal geometry (not necessarily quasi-one- 
dimensional). Subsequent to experiments by Liefrink et 
al ||6[, which demonstrated shot- noise suppression close 
to 1/3 even for conductors much longer than the de- 
phasing length, de Jong and Beenakker [[7| provided 
a semi-classical discussion which showed that the 1/3- 
suppression is insensitive to dephasing. More recent ex- 
periments by Steinbach et al B demonstrated the tran- 
sition from the 1/3-suppression regime in wires short 
compared to an inelastic length through an interaction- 
dominated regime ^j^^ to a regime where shot noise 



is suppressed by inelastic scattering jG,0,n2J. A macro- 
scopic metal exhibits no shot noise [|l3 |. 

Here we investigate the shot noise in mesoscopic diffu- 
sive conductors in a multi-terminal geonretry. Primarily, 
we focus on the interference experiment ||l4], analogous 
to the experiment of Hanbury Brown and Twiss in op- 
tics 1^. Namely, we consider conductor, connected to 
four reservoirs a, f3, 7, and 5 at equilibrium (Fig. 1), and 
discuss three types of experiments. In the experiment A 
current is incident from the probe /3, i.e. ^a = fJ"y = 1^5', 
fJ'f3 — fJ'a = gV ^ fj,\ being the chemical potential of elec- 
trons in the reservoir A. In the experiment B current is 
incident from the probe S: /ia = /i/3 = /i-,; fJ-s ^ ^-a — eV . 
Finally, in the experiment C current is incident from both 
probes (3 and 5: /ia = /i^; /J./3 = /J.^; ^/3 — /ia = eV . The 
current correlation in probes a and 7 is measured in all 
the experiments, Sj{t) = -{AIa{t)AI^{0)), j = A,B,C. 

The general analysis of Ref. Q allows to express these 
quantities in terms of scattering matrices s^"^ , with in- 
dices A and i^ labeling the probes. Thus, for zero fre- 
quency and temperature Em one obtains 



(1) 




with quantities ^^ defined as follows. 

Si = Tr(s+"'3s"/5s+7/357/3). 

52 = Tr(.s+"''s"'^.s+'"'s''''); 

53 = Tr(.s+"'3s"'^s+'^'^s'''^); 

54 = Tr(s+"''s"'^^s+'^^s'''^), 



(2) 



the scattering matrices are evaluated at the Fermi sur- 
face, and the trace is taken with respect to channel in- 
dices. 

Thus, Sc T^ Sa + Sb- experiments A and B are not 
additive due to the interference ternrs S3 and S4. It was 
shown in Ref. pdj that these terms have different signs for 
fermions and bosons; hence we will call them exchange 
terms. We define an exchange contribution as 

AS* = Sc — Sa — Sb- 

It follows from the unitarity of matrices s'*''' that the 
quantities Si and S2, which represent the classical result, 
are positively defined ||l^. At the same time, traces S3 
and S4 can have either sign. This means that exchange- 
interference nray either suppress or enhance the classical 
value. 



In a disordered system all these quantities should be 
averaged over impurity configuration. Naively, one might 
think that due to the phases contained in the quantities 
S3 and S4 these will average to zero, and thus the av- 
erage of the exchange term (AS*) vanishes (here angular 
brackets are used to indicate the disorder average). Be- 
low we explicitly calculate disorder-averaged correlation 
functions Sj, and demonstrate that it is not the case. 
The average exchange correlator (AS) generally has a 
nonzero value. An analysis of the exchange correlator 
for chaotic cavities, reported elsewhere [Q, leads to a 
similar conclusion. 
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FIG. 1. Four-terminal conductors; the disordered area is 
shaded. 

The paper is organized as follows. First, we inves- 
tigate disorder averages of scattering matrices starting 



from Eqs. (g) and using the Fisher-Lee relation which 
connects scattering matrices and Green's functions. We 
then use diagram techniques developed for disordered 
systems to find the ensemble averages. As a simple check 
of the method developed, we give a novel derivation of 
the 1/3-suppression of the two-terminal shot noise for an 
arbitrary (not necessarily quasi-one-dimensional) geom- 
etry, thus confirming the result by Nazarov ^. Then 
we turn to exchange-interference experiments and con- 
sider the two particular four-terminal geometries, shown 
in Fig. 1. We demonstrate that the geometry of Fig. la 
implies a negative exchange correlation, with the quan- 
tity AS" being of the same order of magnitude as correla- 
tors Sa and Sb themselves. In contrast, the cross geom- 
etry of Fig. lb shows a strong suppression of exchange 
effects, and gives a positive sign of the latter, provided 
the motion through the center of a cross is ballistic. Oth- 
erwise the exchange effect is governed by the scattering 
inside the cross center only. 

In the calculations below we disregard electron- 
electron interaction. The latter is known not to produce 
an essential effect on two-terminal shot noise P,|l^l pro- 
vided the wire is short in comparison with the inelastic 
scattering length. We will show that the origin for this is 
that in the ensemble averaged quantities the effect is local 
and electron trajectories enclosing a large area are sup- 
pressed. This explains why the shot noise is not sensitive 
to dephasing. Hence, we believe that electron-electron 
interactions are not important for the exchange effects 
in shot noise. Note, however, that non-linear noise is 
affected by interactions, as was shown recently p9i. In- 
teractions are also expected to affect the frequency de- 
pendence of the shot noise power. 



II. GENERAL FORMALISM AND 
TWO-TERMINAL SHOT NOISE 

We consider a disordered two-dimensional system, con- 
nected to reservoirs by ideal leads. Transverse motion of 
electrons in each lead is quantized, and we assume that all 
leads are wide, i.e. the number of transverse channels at 
the Fermi surface in the lead A is large, N\ — pfW\ 3> 1. 
Here pp is the Fermi momentum, while Wx is the width 
of the lead. 

General relations pc|~[2^] allow one to express scatter- 
ing matrices for an arbitrary geometry through retarded 
and advanced Green's functions of the system. The stan- 
dard procedure [g2[ is as follows. One chooses arbitrary 
cross-sections of the leads C\, and introduces local co- 
ordinates related to these cross-sections (Fig. 2). Since 
nothing depends on the choice of these cross-sections, it 
is convenient to choose them as boundary between disor- 
dered region and leads. One obtains 
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FIG. 2. Contact of a disordered region (shaded) with an 
ideal lead A. 
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Anj2^ M/2 / '^y^ / dy^^G^irx^r^) 

AM^{v„^v,X' Jcx Jc^ 

X {Dxn\){D^fiy)eyi^{-ik„iXx~ iknXy) 

X Xm{y\)xn{yu) (3) 



and 

X Xm(yA)xn(y,.)- (4) 

Here Vm = km/M , M being the effective electron mass. 
The longitudinal wavevectors in the lead A are 



km ^ pp - {-Km/WxY 



1/2 



and those in lead v are denoted by fc„. Furthermore, ri\ 
is the unit vector in the direction a; a, while Xm and Xn 
are wavefunctions of transverse motion in the leads A and 
V respectively; for simplicity we choose them to be real. 
Finally, D denotes a double-sided derivative, 

JDg = JVg - gVf. 

In principle, Eqs. @) and (Q) allow one to average 
arbitrary combinations of scattering matrices over disor- 
der, using the standard diagram technique [g^. It seems 
that the approach outlined here has not so far been used 
for the (analytical) calculation of any physical properties. 
However, it is rather close to the Hamiltonian approach, 
employed extensively for the calculation of conductance 
and conductance fluctuations [p4|-p8[. Below we demon- 
strate that our formalism reproduces the 1/3-suppression 
of two-terminal shot noise; in particular, as a simplest 
check, we also reproduce the Drude formula for conduc- 
tance. 

The rest of the Section is devoted to the two-terminal 
geometry — a diffusive wire of the length L and width W ^ 
connected to two ideal leads a and (3; L,W ^ I, with I 
being the mean free path. For a moment we assume also 
L ^ W, a restriction, which eventually will be lifted. 
We introduce an axis x directed along the wire, < a; < 
L, and an axis y directed across the wire. The general 
expressions (0), (Eh can be rewritten as 



SmiiE) - 



x[-d^,+tkm][-d^,~tkn]G'l{ri,r2)Qz 






X2 — L 



and 






dyiXniyi) I dy2 



X Xrn{y2)[-dxi ~ikn\[-dx2 + ik„ 



2:1=0, 



xG^(r2,n)i:n' 



(6) 



Two-terminal shot noise power 5* = S{lj — 0) can be 
conveniently expressed through scattering matrices eval- 
uated at the Fermi level |E3,P0i, 



S = — ey(Tr [s+'^t's'^P] - Tr [s+"^s 



+a/3 „Q/3 „+a/3 „Q/3] 



27r 



)• (7) 



Note that the first trace on the rhs is related to the con- 
ductance. 



G=|-(Tr[s+"''s"^]}. 
It is convenient to calculate both traces separately. 



A. Evaluation of (Tr [s+"''s"''] >. 

Using Eqs.(ll) and (^), we find for the conductance 



•Ms*'''"^)'j^,^T^„ 



X / dy2dy3Xn{y2)Xn{yz) I dyidy4,Xm{yi)Xm{yi) 

JCc JCfi 

X [ikjyi — Uxi\ [~*fcn ~ tlx2\ [ikn — Oxsl [~*fcm ~ C)xi\ 

x{G^{r,,r2)G^{r,,r,))\:iZ:iZ^^ , (8) 

where the Green's functions are taken at the Fermi en- 
ergy. Since the averaged Green's functions decay on 
scales of the mean free path, the average product of two 
Green's functions, each of them taken in remote points, 
is only due to the diffusion (see e.g. pM): 



(G^(ri,r2)G^(r3,r4)) = / dradrb{G^{ri,ra)) 

X (G^(n, r^)) (G«(r3, n)) {G^'ira, r4))P(r„ n). (9) 

The diffusion propagator P{r,r') is a solution of the 
equation 

- DVlP{r, r') = {2TTVT'^)-'^5{r - r') (10) 

with appropriate boundary conditions (P = at the con- 
tact to the ideal leads; nVP = at the walls). Here 
J/ = M/2'K, D = vfI/2 and r are density of states, diffu- 
sion coefficient, and elastic lifetime, respectively. Under 



the assumption L ^ W the diffusion can be considered 
to be one-dimensional, and the diffusion propagator does 
not depend on y, 

Pix,x') = iMr^DWL)-'ilil'fJ']^^^ H""^,. (ff) 



2,n 



3,n 



R 



R 



l,in 



4,in 



FIG. 3. Diagram for the conductance. The double dashed 
line represents the diffusion propagator. The position and 
the transverse channel number of the points on the surfaces 
Ca {x — 0) and Cp {x — L) are shown. For example, the 
transverse wavefunction Xm is taken at the point j/i. 

Now we insert Eq. (^ into Eq. (g) . The diagram for g 
is shown in Fig. 3. One can approximate the short-ranged 
Green's functions as follows. 



(G^(r,r')) = -— exp 

PF 



IPF 



21 



x5{y~v'). (12) 

Then, integrating over transverse coordinates, we obtain 



M 



IQDt^LW 



2-^ h. 



h 
PF 



dxn.dxh 



X exp[-a::o/;] exp[-(L - Xh)/l]xa{L - Xb)- (13) 

Taking into account that 



\^J_fl + t^ 



Kn 



PF 



2VF, 



we obtain 



9 = ^PfW. 



(14) 



Multiplied by e^/27r, Eq. ( p^ gives the Drude formula, 
as it should be. 



B. Evaluation of (Tr [s+^'^s^^s+'^'^s"'']}. 

The trace of a product of four scattering matrices can 
be written as 

t = (Ti [s+"'^s"^s+"'3s"/3]\ = \ V ^ 



klran 



ku kl km, k' 



k^'l'^m "-n 



X / dy2dy3dyedyjXiiy2)xiiy3)Xniy6)Xniy7) 

J Co, 



X / dyidy^dy<idys,Xk{yi)Xm{yi)x-m[yb)Xk{y%) 

ICp 

X [ikk ~ 9a;i] [-iki - dx-A [ih - dx^,] [-ikm ~ d^^] 
X [ik,n - dxrj [-ikn - dxe] [ikn - dxr] [ikk - dxg] 
x(G^(ri,r2)G«(r3,r4)G^(r5,r6) 



XG- [r7,rs))\x2^x3=xe=x-r=0 



(15) 



Employing Eq. (pi) again, we find the diagrams shown 
in Fig. 4. We omitted all diagrams containing a single 
electron line connecting two different leads, since these 
are exponentially small; the diagrams (a) and (e) contain 
also counterparts, similar to (c) and (d). 



2,1 
3,1 
6,n 
7,n 



I I 
I I 



"TT" 
I I 
I I 



a) 



b) 



X 



c) 



^ 



d) 



1 1 
1 1 



R 



R 



"TT" 
I I 
I I 



l,k 

4,in 

5,in 
8,k 



I I 
I I 



"TT" 
I I 



I I 
I I 



e) 

FIG. 4. Diagrams for the quantity t in the same notations 
as in Fig. 3. The single dashed hne with a cross represents 
impurity scattering. 

The diagrams (b), (c) and (d) turn out to give the 
leading contribution, whereas others carry small fac- 
tors. Thus, for diagram (a) points yi and 2/4 should 
lie not further apart than a mean free path, which due 
to the orthogonality of transverse wavefunctions implies 
k = m. Therefore, the contribution of this diagram is 
suppressed by a factor (ppW)~^ ^ 1. The diagram (e), 
which is topologically equivalent to (b), is suppressed as 
{pfW)~^. Taking into account the explicit form (|ll| ) 
for the diffusion propagator, and integrating over coor- 
dinates Hi and over one of two pair of coordinates in the 
diffusion propagators (those lying close to one of the ends 
of the wire), we arrive at the expression 



t 



/8 



2~1 

V" 1— I 1 + ^^ I / dradrbdrc 
^ k^\ PfJ J 

(16) 



2{ADt^WL)^ 
X drd{L - Xa){L - Xc)xbXdF{ra, Vb, r^ Vd) 



Here F is the Hikami box |^. It is short-ranged (all 
points Va, rb, re, and r^ should be close to each other), 
and in the Fourier-space has the form 

F{qa,qb, qc, qd) = -MT^vj,{2TT)^S{qa + qb + qc + qd) 

X [2{qaqc + qbqd) + {qa + qc){qb + qd)] ■ (17) 

Integration of the Hikami box over the cross-section of 
the wire yields 

/ F{ra,rb,rc,rd)dyadybdycdyd = Mr^VpW [2dx:,d:c, 
+2dj;^dxa + dx^dxb + dx^dxa + dx„dx^ + dx,dxj\ 

x5{Xa-Xb)S{Xa-Xc)S{Xa- Xd). (18) 



Inserting Eq. (|18| ) into Eq. (|lq ) and performing the 
remaining integrations, we obtain 



t = ^PfW = 2g/3 



(19) 



which immediately gives the 1/3-shot noise suppression. 



C. Universality 

Now we lift the requirement L ^ W, but still consider 
a diffusive system, W,L ^ I. The result (|IJ) for g ~ 
(Tr[s+"''s"^]) is equivalent, in fact, to the Drude formula, 
and is therefore valid for an arbitrary relation between W 
and L. In the derivation of i = (Tr[s+"^s"''s+"^s"'']) we 
should now take into account that the diffusion is not one- 
dimensional any more, and write the diffusion propagator 
in the form 



Piry) 



1 



1 






0q(^)0g(r-'), 



(20) 



instead of Eq. (|ll|). Here 4>q(r) and — g are eigenfunc- 
tions and eigenvalues of the Laplace operator with appro- 
priate boundary conditions. In our particular geometry 
one obtains 

with integers n^; > and n^ > 0. It is easy to see 
that the integration over y\ and y% in the diagrams of 
Fig. 4 (b),(c),(d) places a constraint on the wavevec- 
tor n\y of the diffusion propagator connecting these two 
points, n\y — 2k (unless niy = 0). In the same way, the 
other integrations over yi imply other constraints, which 
due to the (5-function in the expression for the Hikami box 
(O) yield a constraint on the channel indices k,l,m,n. 
Therefore all terms with non-zero transverse harmonics 
are small as {pfW)~^. Up to terms proportional to this 
small parameter the result ( |l9| ) is exact. Thus, the 1/3- 
shot noise suppression is, indeed, universal, and does not 
depend on the ratio W/L, provided the system is diffu- 
sive, in accordance with the conclusion of Ref. B. 

To conclude this section, we compare the method used 
above with other derivations of the 1/3-shot noise sup- 
pression 0-0]. As is well known, there exist two princi- 
pally different methods of calculating conductance. One 
can first evaluate conductivity (which is a local quan- 
tity), starting from the Kubo formula, and then, after 
integration over a cross-section one obtains the conduc- 
tance. Alternatively, one can calculate conductance di- 
rectly, starting from the Landauer formula. (In fact, our 
derivation of the quantity g given above is of this kind). 
Both derivations are equivalent, although at intermediate 
stages they have not much in common. 

A similar situation happens in the calculation of shot 
noise. On one hand, one can calculate the microscopic 
correlator of currents, and upon integration over a cross- 
section obtains the shot noise power. The derivation of 
Altshuler, Levitov and Yakovets Q is exactly of this type 
p3[ . It can be generalized to an arbitrary geometry, 
and, in principle, can be used for a broad class of prob- 
lems. The local current correlator contains more infor- 
mation than is necessary for the calculation of the shot 
noise power. The method of Nagaev ||^ and de Jong 
and Beenakker |0], who employ the quantum Langevin 



equation, is somewhat similar, although the equivalence 
between these two approaches is not evident. The gen- 
eralization of the latter approach for a multi-terminal 
geometry does not seem to be quite obvious. 

The derivation of Beenakker and one of the authors 
0, as well as the present method, belong to another, 
scattering (or Landauer) type of approaches. Ref. g] 
derives the shot-noise power with the use of the distri- 
bution of transmission eigenvalues of diffusive wire. This 
proof seems to be the most elegant. However, one should 
not forget that the distribution of transmission eigenval- 
ues itself is derived by sophisticated methods such as the 
DMPK equation Q. Although Nazarov ||] succeeded 
in extending this derivation to the case of an arbitrary 
two-terminal geometry, most probably it can not be gen- 
eralized to multi-terminal case: for conductors with four 
(or more) probes the shot noise is not expressed through 
eigenvalues of the scattering matrix s'^s. The derivation 
given in this paper is more general, and self-contained; it 
does not require the distribution of transmission eigen- 
values, and hence allows a generalization to an arbitrary 
scattering geometry. We repeat that it differs from the 
microscopic calculation of ref. Q. Certainly, although 
these two approaches can not be easily compared at in- 
termediate stages, in the end they should lead to identical 
results for arbitrary scattering geometries. 



III. MULTI-TERMINAL SHOT NOISE 

Now we turn to the exchange-interference experiment 
described in the Introduction. We consider a four- 
terminal geometry (examples are shown in Fig. 1; for 
convenience, we still use the coordinates of Fig. 2), and 
calculate the current correlation for experiments A, B, 
and C. 

The quantity Si = TT{s+°'^s"f^s+'^f^s''P) is determined 
by diagrams of Fig. 4, where now points yi,y4,y5,y8 
belong to the contact with lead /3; points ^2,2/3 and 
2/6,2/7 belong to the contacts with leads a and 7, re- 
spectively. Therefore the diagram (e) is exponentially 
small, while the diagram (a) is suppressed in the param- 
eter {pfWi3)~^. Hence, the quantity Si is given by the 
same diagrams (b),(c), and (d), as the quantity t. Us- 
ing Eq. (|l3) and integrating over the cross-section of the 
leads, we obtain (Fig. 5) 

2^^^ fcl5;;'n ^*'-^'^"^" V pf) \ PF, 

/•O 
,2/„. \,.2 /„, \ / 7 , ,,2 , 



xP[rc, rd)P{re, rf)P{rg, rh)F{rb, r-g, r^, rg 



(21) 



xXfe(ya)Xm(yc) / dynXniVh) I dXadXcdxfdXh 

J C-y J —OQ 

X exp((a;a +Xc + Xf + Xh)/l) dVa-.. drhP{ra, rt) 



Here the points ra,rc,rf, ru are given in the coordinates 
of the contacts (3, (3, a, 7, respectively. 
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FIG. 5. Typical diagram for the quantity Hi. 
In the same way, for the quantity S3 one obtains 



1 fM 






"' 2 V 4 if- kkkk^kn {^^ pf) V^ pf 

klran 

xii + ^)'(^i + ^y|^ dyfxf{yf)J^ dyaxUva) 

X / dyhXniyh) / dycXmiVc) / dxadxcdxjdxh 

JC-, JCs J -00 



X exp((a;a +Xc + Xf + Xh)/l) j dra . ■ ■ drhP{ra, rb) 
xP{r,, rd)Pir,,rf)P{rg, rh)F{rb, r^.n, r,,), (22) 

and the points ra, r^^ r f, r^ are given in the coordinates 
of the contacts /3, (5, a, 7, respectively. Expressions for 
the quantities S2 and S4 can be obtained from Eqs. ( |2l| ) 
and (p3), respectively, by interchanging /? <-> ^. 

Expressions (|2^) and (22) are valid for an arbitrary 
four-terminal geometry and can be used for numerical 
calculations. It is important that not only traces Si and 
S2, as one could expect, but also quantities S3 and S4 
are phase insensitive. Indeed, the electron motion which 
Eqs. ( pi|) and (|2^ ) imply is just the diffusion between 
different leads. No closed paths are formed, except for 
ballistic motion due to the scattering described by the 
Hikami box somewhere in the middle of the sample. Since 
the size of this loop is very small, of order of the mean 
free path, dephasing is not expected to have an effect 
on the exchange noise. Certainly, some effects similar to 
weak localization exist, however, as for conductance [ pl| , 
they are relatively weak (as (pi?/)~^) in comparison with 
the main effect. We do not discuss these effects here. 

To make further progress we have to solve the diffusion 
equation in a given geometry with appropriate boundary 
conditions. We turn now to the two different geometries, 
shown in Fig. 1. 



A. Box geometry 

First, we consider the geometry of Fig. la. We assume 
all leads to be wide, Wx :s> L Then points r^, r^ Vf, 
and r/i are typically far from the lead's boundaries. This 
means that, for example, in the integral over yf one can 
replace the diffusion propagator, P{re,rf) by another 
function P{re,rf), which is also a solution to the dif- 
fusion equation, but with another boundary conditions, 
appropriate for an open surface, 

P(r,r')U=o = 0. 

We do not need to specify boundary conditions for 
P{re,rf) on the other boundaries, since the point r^ 
is typically in the middle of the sample. Consequently, 
we may substitute for all "true" diffusion propagators P 
the functions P, the solution with P = everywhere on 
the boundary, as is appropriate for an open system. The 
solution P is 



P{r,r') = 



E 



1 






X sin — - — sin — sin — — sin — . 

Ijrp Ijrv. J_J^, J_J,f, 



(23) 



Furthermore, the functions P vary considerably on the 
scale of the size of a sample, L^ and Ly. If we assume 
W\ <^ Lx, Ly, the function P{re, rj) in the integral over 
yf may be taken independent of yf. Thus, we obtain 



1 /M 

2 V "2 
L. 



■=1 = :t ( -;r I WaWpW^ / dyadycdxfdxh 



X exp 



y 



Va 



Ly - yc 



I 



Lx ^h 



I III 

xdrddredrgP[Xf3, ya, rb]P[Xp, yc; rd]P[re; Xf,Ya 
xP[rg; Xh, Ky]F(rb, r^, r^, rg) 



dvh 



(24) 



and 



"3 



M 



X exp 



2 

Ly — Va 



WaWpW^Ws / dyadycdxfdxh 



Vc 



"Hi 
I 



Xh 



drh 



III I 

xdrddredrgP[Xp, ya, rb]P[Xs, yc] rd]P[re; Xf,Ya] 
xP[rg; Xh, Y^]F{rb, r^, ra, Vg). (25) 

Here Y^, Xfj, Y^, and X^ denote the positions of the 
corresponding leads. 



We see already from Eqs. ( p4|) and (25) that the re- 
sults are not-universal in the sense that they depend on 
the geometry of the sample. Indeed, within the approxi- 
mation in which we replace P by P, the quantity Si does 
not contain any information on the location and width of 



lead 5; at the same time, it depends essentially on the lo- 
cation and width of other leads. The quantity S2 contains 
information of all leads except /3, whereas both S3 and 
S4 are governed by the geometry of all leads. Therefore 
all ratios S,;/Sj depend essentially on the geometry of 
the sample. This is in contrast with the case of a chaotic 
cavity pq ], where one obtains Si = S2 = — 3S3 = — 3S4 
irrespectively of geometry, provided the leads are wide 
enough. 

Performing the integration and taking into account 
that the remaining sums are converging rapidly for L^ ^ 
Ly (the case we assume from now on), one obtains cum- 
bersome expressions for the quantities S^. In the sym- 
metric case, Lx = Ly = L, W\ = W , Ya = Y^ = Xp = 
Xs = L/2 they simplify. We obtain 



-1 — =^2 
13 — S4 



m 



^/3|^^Ht 



(26) 



with positive constants 

rji = — ^ — 3 — (cosliTT — l)(27rcosh7r — sinhyr) « 0.21 
2 sinh TT 



and 



1 



m 



sinh' TT 



-(27rcosh7r — sinhyr) « 0.03 



It is seen that the exchange effect exists, and has a neg- 
ative sign (i.e. exchange suppresses the result of exper- 
iment C in comparison with the sum of the results of 
experiments A and B). Although the relative value of 
the effect is S3/S1 ^ 0.1, the effect should be clearly 
observable. 



B. Cross geometry 

We consider now the cross geometry of Fig. lb. We 
assume that all arms of the cross have equal P5|] lengths 
L and widths W . For L ':^ W we can consider diffusion 
as one-dimensional. We also assume that the center of 
the cross is described by a reflection coefficient R and a 
transmission coefficient T = (1 — i?)/3 between any two 
different arms. 

The diffusion propagator is a solution of Eq. (|l^) . We 
move to the coordinate system of Fig. lb and fix the point 
r near the origin of the lead a, x c^ L. We introduce 

Pa\{x, x') — P{x, x'), if x' lies in the arm A, 

which is proportional to the time-integrated probability 
of diffusion from point x in the arm a to point x' in the 
arm A. The solution satisfying the boundary conditions 
and the condition of current conservation in the cross. 



^dx'Pa\{x,x')\x'=Q = 0, 



-* aa v-^i ^a.) 



1 {L-x)(Le + 3xa) 



MDt^W 



3 + e 






i'aA(a;,a;A) = ^ — ^^ ^^^ —, X ^ a. (27) 

The constant e, defined as the ratio of diffusion probabil- 
ities, 



^" Pa/3(a;,0)' 

is calculated in the Appendix. The result is 

' 1 + 1{LT)-\1-2T), T:^l/L 
^~ ^ 1{LT)-^, T^l/L 



(28) 



(29) 



Now we substitute Eq. ( p7| ) into the general expres- 
sions (^IJ) and (^2|). Since the area of the cross is negli- 
gible in comparison with the areas of the arms, we can 
neglect the possibility of finding the Hikami box inside 
the cross, and allow it to be situated only in one of the 
arms. Upon integration we obtain 

„ „ I ^^. 3(l + e^) + 4 

^,=^, = -Wp^ (3 + e)4 ' 



=^3 — '^4 



4? e-1 

L^^^(3T^- 



(30) 



Thus, in the case T 3> l|i^^ when the overall transmis- 
sion through the sample is governed by the diffusive arms 
rather than by the center of the cross, one has e ^ 1. The 
quantities Si and S2 are regular for e = 1, and therefore 
assume the finite value, 2i = S2 = (5/192)(pfW^V^)- 
At the same time, the exchange terms S3 and S4 are 
strongly suppressed in the parameter Z/L, S3 = S4 = 
{\IM)(j)fWI'^ I L'^T)(\ - 2T). In the less realistic case 
T <g; l/L (the transmission is determined by the cen- 
ter of the cross) one obtains e » 1. All quantities Si 
are small, since now all channels are nearly closed (cf. 
the situation for two-terminal shot noise [ p9| , p0| ), how- 
ever exchange terms are additionally suppressed in the 
parameter e~^. 

Thus, in the cross geometry of Fig. lb the exchange 
noise (AS") is suppressed in comparison with the regular 
terms {Sa + Sb) irrespectively of the transmission prop- 
erties of the center of the cross. It is also quite remarkable 
that for the cross geometry the exchange contribution is 
positive, although small: the total effect is enhanced by 
the exchange. 



IV. CONCLUSIONS 

We have investigated shot noise in diffusive conductors 
on the basis of Eq. (||) and the Fisher-Lee relation, which 
expresses scattering matrices through advanced and re- 
tarded Green's functions. In this way, one can reduce 



disorder averages of various combinations of scattering 
matrices to standard diagram technique for Green's func- 
tions [^ . Although this approach resembles previously 
published calculations of conductance and conductance 
fluctuations |24-E3], we believe it to be more transparent. 
We are not aware of any applications of this approach to 
noise problems. 

As a check of the method, we first reproduced the 
1/3-shot noise suppression in the two-terminal geome- 
try and confirmed the statement of Ref. Q that it is in 
fact super-universal and holds for an arbitrary relation 
between length and width of a wire, provided the system 
is diffusive. Our proof bears some similarity with other 
ones existing in the literature |0-B; however, it is novel, 
and a direct equivalence to any of the existing proofs is 
not evident (see the discussion in the end of Section 2). 

Then we turned to the multi-terminal geometry and 
investigated the interference experiment, similar to the 
Hanbury Brown and Twiss experiment known in optics 
|15[ . We obtained general expressions for scattering ma- 
trix combinations (p^), (p3), determining noise intensi- 
ties (IT]); then we investigated them for the two different 
geometries of Fig. 1. 



a 



-<>' 



FIG. 6. Typical electron trajectories, contributing to the 
quantity E3. Solid lines denote ballistic propagation (de- 
scribed by averaged single-particle Green's function), and 
dashed lines denote diffusive propagation (described by the 
diffuson P). 

The important point we make is that the exchange 
effect, even when averaged over disorder, does not van- 
ish. The reason is that typical electron trajectories, con- 
tributing to all averaged traces of scattering matrices, 
considered above (i.e. quantities g and t for the two- 
terminal geometry, and S.^ in the four-terminal case) do 
not contain large closed loops. In particular, it is valid 
for the "exchange" traces S3 and S4. A typical trajec- 
tory for the quantity S3 is shown in Fig. 6. It is a di- 
rect translation of diagrams contributing to this quan- 



tity. The electron motion is essentially diffusion between 
different leads with ballistic propagation (described by 
disorder-averaged single-particle Green's function) close 
to the leads and somewhere in the middle of the sam- 
ple (the later motion described by the Hikami box in 
Eq. (P^)). Thus, closed loops are related to ballistic mo- 
tion over distances of an elastic scattering length only, 
and therefore neither the shot noise in two-terminal con- 
ductors nor the shot noise in multi-terminal structures 
should be sensitive to dephasing. 

Another observation is that exchange corrections are 
not universal, in contrast to what is found in the chaotic 
case Q: the ratio {AS)/{Sa + Sb) depends on the ge- 
ometry of a sample in an essential way. Even the sign of 
the effect may change: for the box geometry of Fig. la 
it is negative, i.e. interference suppresses the total ef- 
fect, while for the cross geometry (Fig. lb) interference 
enhances the effect (although weakly). 

The results obtained for the cross geometry allow us 
to make predictions for experiments in real systems. In- 
deed, we found that the exchange contribution is sup- 
pressed strongly with respect to the average noise in- 
tensities (Sa) and {Sb}- This result was obtained by as- 
suming that the intermediate scattering, described by the 
Hikami box, does not happen in the center of the cross, 
i.e., strictly speaking, for ballistic propagation through 
the center. In more complicated situations the entire 
exchange effect will be determined by properties of the 
center of the cross. If the motion within the center is dif- 
fusive, one can apply the results obtained above for the 
box geometry. The total exchange effect is expected to 
be negative. However, since the arms of the cross (which 
correspond to disordered leads in the real experiments) 
contribute to the intensities (Sa) and (Sb), but not to 
the exchange contribution, the latter will still be sup- 
pressed, if disorder extends far into leads. Finally, if the 
center of the cross is a chaotic cavity, one may use the 
results of Ref. [Q. The exchange contribution in the 
chaotic cavity separated from ideal leads by high barri- 
ers (disordered arms play the role of these barriers) is 
positive: the interference enhances the effect. 



ACKNOWLEDGMENTS 

We thank S. van Langen, who calculated the exchange- 
interference correlator for chaotic cavities, for useful dis- 
cussions. The work was supported by the Swiss National 
Science Foundation. 



number of scatterers in each arm is A^ = L/l. Each scat- 
terer is described by transmission t — 1/2 and reflection 
r = 1/2 probabilities. We denote the carrier flux densi- 
ties in the arm a between sites n and n + 1 away from the 
center of the cross by a„, and the flux towards the center 
of the cross by 6„. Corresponding amplitudes in other 
arms are denoted by a'^ and 6'j (Fig. 7). The total flux 
at each site is given by p„ = a„ -I- 6„, and p'^ = a'^ + b'^. 
The coefficient e can be expressed as e = Po/ p'o- 



/ • 



^2 ^1 ^0 



• • 



^2 ^1 ^0 



\ bj b2 



• • 



^0 ^\ ^2 



FIG. 7. Discrete diffusion model. 

The diffusion equation implies that all densities should 
be linear functions of n; furthermore, matching condi- 
tions at each scatterer require 6n-i — a,n', b'^-i = o,'^- 
Thus, we write 

an^A + B{n~l), a'^ = A' + B'{n - 1), (Al) 



bn^A + Bn, b' =A' + B'n. 






(A2) 



The four constants A, B, A' , B' obey four equations: 

1) Boundary condition for the arm f3: b'j^ — 0. 

2) Matching conditions at the center of the cross: 



and 



a'o = Tba + 2Tb'„ + Rbo, 



aa = 3T&Q + Rbo. 



(A3) 



(A4) 



The fourth equation is Eq.dTfl), however, it is not re- 
quired for the calculation of the constant e. We obtain 



2(Ar + T-i) -3 
2NTi ' 



(A5) 



and the limiting cases given by Eq. (|29|) follow immedi- 
ately. 



APPENDIX 

To find the coefficient e defined by Eq. (|28| ) it is in- 
structive to consider a discrete model of diffusion Q. 
Each arm is modeled by a one-dimensional array of scat- 
terers, placed at a distance / from each other; the total 
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